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Mixed Map Labeling 

Maarten Ldffler* Martin Ndllenburg^ Frank Sfaa/s* 


Abstract 

Point feature map labeling is a geometric problem, in which a set of input points must be labeled with a set of 
disjoint rectangles (the bounding boxes of the label texts). Typically, labeling models either use internal labels, 
which must touch their feature point, or external (boundary) labels, which are placed on one of the four sides of the 
input points’ bounding box and which are connected to their feamre points by crossing-free leader lines. In this 
paper we study polynomial-time algorithms for maximizing the number of internal labels in a mixed labeling model 
that combines internal and external labels. The model requires that all leaders are parallel to a given orientation 
9 G [0, 2tt), whose value influences the geometric properties and hence the running times of our algorithms. 


1 Introduction 

Annotating features of interest in information graphics with textual labels or icons is an important task in 
information visualization. One classical application, whose principles easily generalize to the labeling 
of other illustrations, is map labeling, where labels are mostly placed internally in the map. Common 
cartographic placement guidelines demand that each label is placed in the immediate neighborhood of its 
feature and that the association between labels and features is unambiguous, while no two labels may 
overlap each other [13,21]. Point feature labeling has been studied extensively in the computational 
geometry literature, but also in the application areas. It is known that maximizing the number of 
non-overlapping labels for a given set of input points is NP-hard, even for very restricted labeling 
models [9,18]. In terms of labeling algorithms, several approximations, polynomial-time approximation 
schemes (PTAS), and exact approaches are known [1,7,16,22], as well as many practically effective 
heuristics, see the bibliography of Wolff and Strijk [23]. If, however, feature points lie too dense in the 
map or if their labels are relatively large, often only small fractions of the features obtain a label, even in 
an optimal solution. 

An alternative labeling approach using external instead of internal labels is known as boundary 
labeling in the literature. This labeling style is frequently used when annotating anatomical drawings and 
technical illustrations, where different and often small parts are identified using labels and descriptive texts 
outside the actual picture, which are connected to their features using leader lines. While the association 
between points and external labels is often more difficult to see, the big advantages of boundary labeling 
are that even dense feature sets can be labeled and that larger labels can be accommodated on the margins 
of the illustration. Many efficient boundary labeling algorithms are known. They can be classified by fhe 
leader shapes that are used and by the sides of the picture’s bounding box that are used for placing the 
labels [2,5,6,10,12,15,20]. 

The combination of internal and external labeling models using internal labels where possible 
and external labels where necessary seems natural and has been proposed as an open problem by 
Kaufmann [14]; however, only few results are known in such hybrid or mixed settings. In a mixed 
labeling, the final image will be an overlay of the original map or drawing, a collection of internal text 
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{a) 9 = TT 


(b) 0 = 0 


(c) 6 — 7r/3 


Fig. 1: A sample point set with mixed labelings of three different slopes. In (a) five external labels 
are necessary, whereas (b) and (c) require only four external labels. The slope in (c) yields 
aesthetically pleasing results. 

labels, and a collection of leaders leading out of the image. Depending on the application, we may wish 
to forbid intersections between some or all of these layers. When no additional intersection constraints 
are imposed, the problem reduces to classical internal map labeling. Under the very natural restriction 
that leaders cannot intersect internal labels, the internal labels need to be placed carefully, as not every 
set of disjoint internal labels creates sufficient gaps for routing leaders of the prescribed shape from all 
remaining feature points to the image boundary. Ldffler and Nollenburg [17] studied a restricted case 
of hybrid labeling, where a partition of the feature points into points with internal fixed position labels 
and points with external labels to be connected by one-bend orthogonal leaders is given as input. They 
presented efficient algorithms and hardness results, depending on three different problem parameters. 
Bekos et al. [3] studied a mixed labeling model with fixed-position internal labels and external labels on 
one or two opposite sides of the bounding box, connected by two-bend orthogonal leaders. Their goal 
is to maximize the number of internally labeled points, while labeling all remaining points externally. 
Polynomial and quasi-polynomial-time algorithms, as well as an approximation algorithm and an ILP 
formulation were presented. 


Contribution. In this paper, we extend the known results on mixed map labeling as follows. We present 
a mixed labeling model, in which each point is assigned either an axis-aligned fixed-position internal 
label (e.g., to the top right of the point) or an external label connected with a leader of slope 6, where 
9 G [0, 2tt) is an input parameter defining the unique leader direction for all external labels, measured 
clockwise from the negative x-axis (see Fig. 1). In this model, we present a new dynamic-programming 
algorithm to maximize the number of internally labeled points for any given slope 9, including the left- 
and right-sided case {9 = 0 or 9 = vr), which was studied by Bekos et al. [3]. While for the right-sided 
case Bekos et al. provided a faster 0{n log^ n)-time algorithm, where n is the number of input points, our 
algorithm improves upon their pseudo-polynomial 0(n^°®”+^)-time algorithm for the left-sided case. We 
solve this problem in 0(n^(log n + 6)) time, where 6 = min{n, l/dmin} is the inverse of the distance 
dmin of the closest pair of points in V and expresses the maximum density of V (Section 2). In the 
general case it turns out that the set of slopes can be partitioned into twelve intervals, in each of which the 
geometric properties of the possible leader-label intersections are similar for all slopes. Depending on 
the particular slope interval, the amount t(n, 6 , 9) of “interference” between sub-problems varies. This 
significantly affects the algorithm’s performance and leads to running times between 0{rfi logn) and 
0(n^(logn -f t(n, d, 9))) = (Section 3). Moreover, we can use our algorithm to optimize the 

number of internal labels over all slopes 9 at an increase in running time by a factor of 0{v?‘), as is shown 
in Section 4.1. 


Problem Statement. We are given a map (or any other illustration) M., which we model for simplicity 
as a convex polygon (this is easy to relax to larger classes of well-behaved domains), and a set V of n 


2 











points in A4 that must be labeled by rectangular labels (the bounding boxes of the label texts). In addition, 
we are given a leader slope 6 G [0, 27r). For simplicity we assume that 6 is none of the slopes debned 
by two points in V. We discuss in Section 4.4 how to remove this resttiction. There are two choices for 
assigning a label to a point p G "P: either we assign an internal label Xp on AI in a one-position model, or 
an external label outside of M. that is connected to p with a leader 7 ^. An internal label Xp is a rectangle 
that is anchored at p by its lower left corner. A leader 7 ^ is a line segment of slope 9 inside M.\ it may 
bend to the horizontal direction outside of A4 in order to connect to its horizontally aligned label, see 
Fig. Ic. So in this model, the labeling is fixed once the choice for an internal or external label has been 
made for each point p G P. For a valid label assignment we require that (i) the internal labels do not 
overlap each other or the leaders, and that (ii) the leaders themselves do not intersect each other. Figure 1 
shows valid mixed labelings for three different slopes. 

Given a set of points P ^ P, let A(P) = {Ap | p G P} denote the set of (candidate) labels 
corresponding to the points in P and let r(P) = { 7 p | p G P} denote the set of (candidate) leaders 
corresponding to the points in P. A labeling of P is a partition of P into sets X and E, the points in X 
labeled internally, the points in £ labeled externally, such that no two labels in A(X) intersect, no two 
leaders in r(X) intersect, and no label from A(X) intersects a leader from r(X). 

For ease of presentation we hrst assume that all labels have the same size, which, without loss of 
generality, we assume to be 1 x 1. Hence, an internal label Xp is a unit square with its bottom left corner 
on p. This may be a realistic model in some settings (e.g., unit-size icons as labels), but generally not all 
labels have the same size. We will sketch how to relax this restriction in Section 4.3. 

Each leader 7 p can be split into an inner part (or inner leader), which is a line segment of slope 9 
from p to the intersection point with the boundary of M., and an outer part (or outer leader) from the 
boundary of Jvi to the actual label. We focus our attention on the inner leaders as they determine how P 
is separated into different subinstances. Hence we can basically think of the leaders as half-lines with 
slope 9. For completeness, we explain a simple method of routing the outer leaders in Section 4.2. 

It is well known that in general not all points in P can be assigned an internal label. The corresponding 
label number maximization problem is NP-hard [9, 18], even if each label has just one candidate 
position [17]. If, however, all labels have the same position (e.g., to the top left of the anchor points) 
and no input point may be covered by any other label, the one-position case can be solved efficiently by 
first discarding all labels containing an anchor point and then applying a simple greedy algorithm on the 
resulting staircase patterns [17]. On the other hand, it is also known that any instance can be labeled with 
external labels using efficient algorithms [4,5]. Mixed labelings combine both label types and sit between 
the two extremes of purely internal and purely external labeling [3, 17]. Here we are interested in the 
internal label number maximization problem, which was first studied for 9 G {0, vr} by Bekos et al. [3]: 
Given a map M., a set of points P in Af and a slope 9 G [0, 27r), we wish to find a valid mixed labeling 
fhaf maximizes fhe number |X| of infernally labeled poinfs. 


2 Leaders from the Left 

We sfarf wifh fhe case fhaf 0 = 0, i.e., all leaders are horizonfal half-lines leading from fhe poinfs fo fhe 
leff of M.. Our approach for maximizing fhe number of infernal labels is fo process fhe poinfs in P from 
right to left and to recursively determine the optimal rightmost unprocessed point p to be assigned an 
external label. Since no leader may cross any internal label, the leader 7 ^ decomposes the current instance 
left of p into two (almost) independent parts, one above 7 ^ and one below. As it turns out, a generic 
subinstance can be dehned by an upper and a lower leader shielding it from the outside and additional 
information about at most one point outside the subinstance. The problem is then solved using dynamic 
programming. 
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(b) 


Fig. 2: The slab S{£, u) in yellow, the region S{i, u) and its points from V in purple, and the region E{i) 
in blue. 


2.1 Geometric Properties 

Let p = {px,Py) he a point in the plane and let Lp = {q \ < Px} and Rp = {q \ qx > Px} denote 

the half-planes containing all points strictly to the left and to the right of p, respectively. Analogously, 
we define the half-planes Tp and Bp above and below p, respectively. Let S{i, u) = TiD Bu denote the 
horizontal slab defined by points £ and u (with £y < Uy), and let S{£, u) = S{£, u) H LiD Lu denote 
the set of points in this slab that lie to the left of both £ and u, see Fig. 2. We define Pi^u as the subset 
of V in S{£, u) including £ and u, i.e., Pi^u = 'Pfl {S{£, u) U {£, u}). With some abuse of notation we 
will sometimes also use Lp, Rp, Tp, and Bp to mean the subset of V that lies in the respective half-plane 
rather than the entire half-plane. 

Recall that 6 = min{n, 1 /dmin} is a parameter that captures the maximum density of V as the inverse 
of the smallest distance dmin between any two points in P. We can use S to bound the number of points 
in a unit square that may be labeled internally. 

Lemma 1. At most 0{5) points in any unit square have a label A that does not contain another point in 

V. 

Proof. Let Ehea unit square and let E'p = E r\V he the input points in E. Since the label Xp of every 
point p G Ep is a unit square anchored by its lower left comer at p, no other point q £ V may lie to 
the top-right of p —otherwise p must be labeled externally. Hence any set of points in Ep whose labels 

do not contain another point of V must form a sequence p^^\p^‘^\ ... such that px^ < px'^ and 
fi) (i) 

Py' > py for any i < j. Since the minimum distance of any two points is dmin we immediately obtain 
that Ep contains at most 0(d) points that can be labeled internally. □ 

Next, we characterize which leaders or labels outside of S{£, u) can interfere with a potential labeling 
of Pi^u assuming £ and u are labeled externally. 

Lemma 2. Let £,u € V, let (I' , £'), with £,u € £' be a labeling ofPf^u- There is no point in T„ U Bi 
whose leader intersects a label from A(X') and there is no point in T^ whose label intersects a label from 
A{I'). 

Proof In any labeling of „ there are no intersections between labels and leaders. In particular, no label 
Xp for p £ T' intersects or 7 ^. It follows that all labels for T' lie inside the slab S{£, u). By definition, 
all leaders for S' also lie inside S{£, u). Since leaders of points in Tu U Bi do not intersect S{£, u), no 
such leader can intersect a label from I'. Moreover, all labels are anchored by their bottom left corner, 
and hence all labels of points in T^ lie above u and do not intersect S {£, u) . Thus, no label for a point in 
can intersect a label in A (X'). □ 

It is not true, however, that labels for points in Bi cannot intersect labels for Pi^u- Still, the influence 
of Bi is very limited as the next lemma shows. Let E{£) denote the open unit square with top-left corner 
£, i.e., E{£) = Rif] BiCi Lr Cl Th, where r = {£x + 1, £y) and b = {£x,£y — 1)- See Fig. 2b. 
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Fig. 3: (a) u, r) expresses the maximum number of points in S{£, u) (purple), that can be labeled 
internally in the depicted situation, (b) The rightmost point p that is labeled with an external label 
decomposes the problem into two subproblems (the orange and blue points). 


Lemma 3. Let i,u € V, let {X' with i,u G £' be a labeling of and let (X",f") denote a 
labeling ofV D BiU {f\ with i G E" . There is at most one point p G X" whose label may intersect a 
label ofX', and p G E(£). 

Proof. Since I ^ 8' r\ 8" we know that no label in A(X') or h.{X") intersects 7 ^. Thus 7 ^ serves as a 
separation line between the labels A(X') and A(X"). Let V" C X" denote the set of points whose labels 
intersect a label of X', and let A" = A('P") denote the corresponding set of labels. We first argue that 
V" C E{i). Then we argue that V" can contain at most one point. 

The labels in A" do not intersect 7 ^, hence they lie strictly right of i. Thus, V" C R^. All points in X' 
lie to the left of £, and their labels have width one. All labels in A" intersect such a label, thus all points 
in X" must lie in L^, where r = {£3; + l,£y). All labels in A" intersect a label from a point in X'. Thus, 
all labels in A" intersect the horizontal line containing 7 ^. Since all labels have height one, it follows that 
all points in V" lie in Th, where b = {£x,£y — 1)- By definition the points in V" lie in Bi. So we have 
V" X R^nLrnTbnBi = E{£). 

The labels in A" are pairwise disjoint and all intersect the top side s of E{£). Since the length of s is 
smaller than one, each label in A" has width exactly one, and all labels lie in R^ it follows that there can 
be at most one label in A". Thus, there is also at most one point in V”. □ 

From Lemma 2 and Lemma 3 it follows that if £ and u are labeled externally, there is at most one 
point r below £ that can influence the labeling of the points in S{£, u). 

2.2 Computing an Optimal Labeling 

We define <h(£, u, r), with £,u £ V, and r G E{£) U {X} as the maximum number of points in S{£, u) 
that can be labeled internally, given that 

(i) the points £ and u are labeled externally, 

(ii) all remaining points in S{£, u) \ S{£, u) have been labeled internally, and 

(Hi) point r is labeled internally. If r = X then no point in E{£) is labeled internally. 

See Fig. 3(a) for an illustration. Furthermore, given £, r, and a point p G Tiwe define g{p, £, r) to be 
the topmost point in E{p) n (T^ U {r}) if such a point exists. Otherwise we define p(p, £, r) = X. 

Lemma 4. For any £,u £ V, and r £ E{£) U {X}, we have that <!>(£, u, r) = \S{£, tt)|, or ^{£, u, r) = 
\Rp n S{£, u)| X r) X u, r'), where p is the rightmost point in S(£, u) with an external label 

and r' = p(p, £, r). 

Proof Let (X*,8*) be an optimal labeling of S{£, u) that satisfies the constraints (i)-(iii) on ^{£, u, r), 
i.e., <l>(f, u, r) = \X*\. In case 8* = 0, we have u, r) = |S'(£, tt)| and the lemma trivially holds. 
Otherwise, there must be a rightmost point p £ 8* with an external label. Consider the partition of 
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Z* at point p into the lower left part B* = Bp n Lp n X*, the upper left part T* = TpH LpCiZ*, and 
the right part R* = RpHZ*, see Fig. 3h. We show that \R*\ = \Rp n S{i,u)\, \B*\ = r), and 

|r* I = <h(p, u, r'), which proves the lemma. 

Since p is the rightmost point with an external label it follows that all points in S{£, u) right of p are 
labeled internally. Hence, R* = RpH S{ 1 , u). 

Next, we observe that Cb = {B *, S(i,p) \ B*) as a sub-labeling of (X*, £*) forms a valid labeling 
of S{i, p), so by Lemma 3 there is at most one point r below i that can influence the labeling of S{i, p). 
This point f, if it exists, lies in E{ 1 ). By constraint (Hi) point r lies in E(i) or r = X and no point in 
E{£) is labeled internally, and thus r can be the only point in E(£) labeled internally, i.e., f = r. So, we 
have that (i) £ and p are labeled externally, (ii) all points in S{£,p) \ S{l,p) are labeled internally, and 
(in) point r is the only internally labeled point in E{£). Thus the definition of <1> applies and we obtain 
|B*| < r). 

Lemmas 2 and 3 together imply that any labeling of S{£,p) is independent from any labeling of 
S{p, u). Thus, it follows that £5 is an optimal labeling of S{£,p) (given the constraints), since otherwise 
(X*,X*) could also be improved. Thus |B*| > ^(£,p,r) and we obtain |B*| = <!>((', p,r). 

Finally, we consider the upper left part T* ■ By Lemma 3 there is at most one point r' in Bp with an 
internal label that can influence the labeling of S{p, u) and we have r' G E(p). We need to show that 
r' = q{p, £, r). Then the rest of the argument is analogous to the argument for B*. 

We claim that r' is the topmost point in E(p) n (T^U {r}). Assume that r' 0 T^, which means r' G B^. 
We know that 7^ does not intersect A^' and hence r' G R^. This means that r' G B(p) H n =: X and 
since p lies to the top-left of I we have C B(£). By definition r is the only point with an internal label 
in E{£) and hence r' = r. So if r' 7^ r we have r' G B(p) H T^. Now assume that r' ^ q{p, £, r). Then 
there is another point q G B(p) H above r' . This point q must be labeled externally since no two points 
in E (p) can be labeled internally. This is a contradiction since by definition p is the rightmost externally 
labeled point in S(£, u) and by constraint (ii) all points in S{£, u) \ S{£, u) are labeled internally. So 
indeed r' = p(p, £, r) and the same arguments as for B* can be used to obtain |r* | = <h(p, u, r'). □ 

Let £, u £ V, and p G S{£,u). We observe that \S{£,p)\ and |5(p, tt)| are strictly smaller than 
|5(£, u)|. Thus, Lemma 4 gives us a proper recursive definition for fi); 

u, r) = max u)), 

max { 4 '(Bp n S{£, u)) + r) -f ^{p, u, q{p, £, r))} }, 

p^S(t^u) 

where 

{ |p| if all labels in A(B U {r} U (S{£, u) \ S{£, u))) are pairwise disjoint, and 
their intersection with and 7 ^ is empty, 

—00 otherwise. 

We can now express the maximum number of points in V that can be labeled internally using 4>. We 
add two dummy points to V that we assume are labeled externally: a point poo that lies sufficiently far 
above and to the right of all points in V, and a point p-00 below and to the right of all points in V. The 
maximum number of points labeled internally is then 4>(p_oo,Poo, X)- 

We use dynamic programming to compute ^(£, u, r) for all u G B U {poo,P-oo} with ly < Uy and 
r £ E(£) U {X}. By finding the maximum in a set of linear size, each value 4>(f, u, r) can be computed 
in 0(n) time, given that the values u' ,r') for all subproblems have already been computed and 
stored in a table and the relevant values for the functions g and / have been precomputed. There are 
0(n) choices for each of £ and w, further there are 0{5) choices for the point r given £ since r is labeled 
internally and we know from Lemma 1 that there are at most 0(8) points in E(£) as candidates for an 
internal label. This results in an 0(n^5) time and 0(n‘^6) space dynamic-programming algorithm. We 
show next that the preprocessing of g and / can be done in 0(n^ log n) time. 


6 





To compute u, r) we actually have to compute q{p, £, r) and r) := "^{Rp H S{£, u)) for 
all points p G S{i,u). We can preprocess all points in "P in 0(n log n) time, such that we can compute 
each g{p, f, r) in 0(1) time as follows. First, we compute and store for each point p £ V the topmost 
point Qp £ V in E{p). This requires n standard priority range queries that take 0{n log n) time in total 
using priority range trees with fractional cascading [8, Chapter 5]. To compute p(p, £, r) we then check 
if qp lies above 1. If it does, we have g{p, £, r) = qp. Otherwise, the only candidate point for g{p, £, r) 
is r and we can check in 0(1) time if r lies in E{p). This takes 0(1) time for each triple {p, i, r) and 
0{'n?5) time in total. 

Next, we fix £ and u, and compute a representation of tk' in O(nlogn) time, such that for each 
p £ S{£, u) and r £ E{£) U {_L} we can obtain 'k'(p, r) in constant time. 

We start by computing the values 'k'(p, -L), for all p. We sweep a vertical line from right to left. That 
is, we sort all points in S{£, u) by decreasing x-coordinate, and process the points in that order. The 
status structure of the sweep line contains the number of points N in S{£, u) right of the sweep line, and a 
(semi-)dynamic data structure T, which stores the labels from the points right of the sweep line, and can 
report all labels intersected by an (axis-parallel) rectangular query window. All labels are unit squares, so 
\r intersects a label \q if and only if contains a corner point of Ag. Furthermore, we only ever insert 
new labels (points) into T, thus it suffices if T supports only insert and query operations. It follows that 
we can implement T using a semi-dynamic range tree using dynamic fractional cascading [19]. In this 
data structure insertions and queries take O(logn) time. 

When we encounter a new point p, p ^ {£, u} we test if the label of p intersects any of the labels 
encountered so far. We can test this using a range query in the tree T. If p £ S {£, u) we also explicitly 
test if Xp intersects -yi or 7„. If there are no points in the query range Ap, and Ap does not intersect 7^ or 
7 u we report -L) = N, increment N (if applicable), and insert the corner points of Ap into T. If the 
query range Ap is not empty, it follows that '^'{p', -L) = — 00 , for p' = pas well as for any point to the 
left of p. Hence, we report that and stop the sweep. Our algorithm runs in 0{n log n) time: sorting all 
points takes 0{n log n) time, and handling each of the 0{n) events takes 0(log n) time. 

Now consider a point r £ E{£). We observe that for all points p right of r, we have that 'k'(p, r) = 
T''(p, _L) = 0 since r is right of all points in S{£, u). Consider the points left of r ordered by decreasing 
x-coordinate. There are two options, depending on whether or not Ar intersects the label Ap of the current 
point p. If \r intersects Ap, we have T''(p, r) = —00 as well as '^'{p' ■, r) = —00 for all points p' left of p. 
If Ar does not intersect Ap we still have T''(p, r) = T''(p, _L). We can test if Ar intersects any other label 
using a range priority query with Ar in (the final version of) the range tree T. We need 0((5) such queries, 
which take 0(log n) time each. This gives a total running time of 0{n log n). We then conclude: 

The above algorithm can also be used when the leaders have a slope 0/0. However, the data 
structure T that we use is fairly complicated. In this specific case where 0 = 0, we can also use a much 
easier data structure, and still get a total running time of 0(n log n). Instead of using the semi-dynamic 
range tree as status structure, we use a simple balanced binary search tree that stores (the end-points of) a 
set of vertical forbidden intervals. When we encounter a new point p, we check if py lies in a forbidden 
interval. If this is the case then Ap intersects another label. Otherwise we can label p internally. This set 
of forbidden intervals is easily maintained in O (log re) time. 

We use this algorithm for every pair {£, re). Hence, after a total of 0{n^ log re) preprocessing time, 
we can answer T'/p, r) queries for any p and r in constant time. This yields the following result, which 
improves the previously best known pseudo-polynomial 0(re^°§’^+^)-time algorithm of Bekos et al. [3] 
for the left-sided case 0 = 0. 


Theorem 5. Given a set V of n points, we can compute a labeling ofV that maximizes the number of 
internal labels for 9 = 0 in 0{n^ logre -f n^5) time and 0(n^d) space, where 6 = min{re, 1/dmmjfor 
the minimum distance dmin in V. 
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Fig. 4: (a) There may be more than one point “below” I with an internal label if the leaders arrive from 
the bottom left, (b) For other directions there is also a region F{u) “above” the subproblem that 
can influence the labeling of S{£, u). 


3 Other Leader Directions 

For other leader slopes 0 7 ^ 0 we use a similar approach as before. We consider a sub-problem u) 
defined by two externally labeled points I and u. We again find the “rightmost” point in the slab labeled 
externally. This gives us two sub-problems, which we solve recursively using dynamic programming. 
However, there are four complications: 


• The region E{1) containing the points “below” the slab S{1, u) that can influence the labeling of 
S{1, u) is no longer a unit square. Depending on the orientation, it can contain more than one point 
with an internal label. See Fig. 4(a). 

• In addition to the region E{£), which contains points that can interfere with a subproblem from 
below, we now also need to consider a second region, which we call E{u), containing points whose 
labels can interfere with a subproblem from above. See Fig. 4(b). 

• The labels of points in S{£, u) are no longer fully contained in the slab S{£, u). Hence, we have to 
check that they do not intersect with leaders of points outside S{£, u). See Fig. 4(b). 

• The regions E {p) and E {p) are no longer strictly to “the right” of p. Hence, for some sub-problems 
we may have already decided (by definition of the sub-problem) that a point q £ E{p) that lies 
“left” of p is labeled internally. Hence, we can no longer choose q to be the rightmost point in 
S{£,p) to be labeled externally. 

We start by explicitly finding the points in V whose internal labels contain other points. We are forced 
to label these points externally. It is easy to find those points in O(n^) time in total. Let Vx denote this 
set of points. Additionally, we spend 0{n'^) time to mark each point if its label intersects r{Vx)- Hence, 
for each point we can determine in constant time if it intersects T{Vx)- For ease of notation we will write 
V to mean the set V \ Vx in the remainder of this section. 

Let 6 be the given orientation for the leaders. Consider conceptually rotating the coordinate system 
such that orientation 9 corresponds to the negative x-axis as in the previous section, and let Bp,Tp, Lp, 
and Rp denote the points in V in the bottom, top, left, and right half-planes bounded by p with respect to 
this coordinate system. Analogously, we define S{£, u) = Tg Cl and S{£, u) = S{£, u) n n L^. 

Our goal is again to bound the number of different labelings of the points in Bi and that can 
influence the labeling of S{£, u). We will show that whether a labeling of Be influences the labeling of 
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S{1, u) depends only on a small subset of the points in Bi. We refer to a labeling of restricted to those 
points as a configuration of B^. The same holds for the points in Tu- 

3.1 Bounding the number of Configurations 

We start by bounding the number of configurations of Bg. Let E{i) denote the bottom infiuence region of 
i. That is, the points in B^ “below” the slab S{£, u) whose label can intersect a label of a point in S{i, u). 
Fig. 5 shows the regions E{1) for various orientations of the leaders. 

Similarly, we can define a region E'{€) C B^ such thaf fhe leaders of poinfs in E'{£) can infersecf a 
label of a poinf in S{£, u). 

Lemma 6. For a subproblem S{i, u) the size of the bottom infiuence region E(£) is at most 1 x e{0) or 
e{6) X 1, where 

0 if 9 G [Svr/d, 57r/4] 

1 G (57r/4,37r/2) 

0 if 9 = 37r/2 
3 if 9 G (37r/2, 77r/4) 

2 G [77r/4,27r). 


Proof We prove fhis by case disfincfion on 9. 
case 0 = 0. See Lemma 3. 

case 9 G (0, 7 r/ 4 ). All labels in E{S{£, u)) lie in Ti and in Lr, where r = {£x + 1, £y). It follows that 
E{£)C Til n L^i, where £' = {£x,£y — 1) and r' = {rx,ry — 1). Furthermore, it is easy to see that 
E{£) c Re, and that E{£) c Be- Hence, E{£) c T^/ Ci ReCi Lr' n Be. 

Since the leaders are sloped downwards it follows that the height of E{£) is at most one. The maximum 
width of E{£) is realized by £ and the intersection point p of the lines bounding Be and Lr'. Using that 
0 G (0, 7 r/ 4 ) basic trigonometry shows that the width is at most two. 
case 0 G [ 7 r/ 4 , 7 r/ 2 ). Similar to the previous case. However, now the width is determined by the 
intersection between and Be. From basic trigonometry it then follows that the width is at most one. 
case 0 = tt/2. Forlabels from £'(^) to intersect A(5 (£,u)) we need £'(^) C T^/, where = {£x,£y — l). 

However, to avoid intersecting je we need E{£) C Bp. It follows that E{£) = 0. 
case 0 G (7r/2, 37r/4). Using similar arguments as before it follows that C BpnLe'FBenLenRe', 
where £' = {lx — l,£y — 1). Since E{£) (Z LeCi Re' the width is at most one. Basic trigonometry 
again shows that the height is at most one. 

case 0 G [ 37 r/ 4 , 57 r/ 4 ]. For the sub-case 0 G [ 37 r/ 4 , tt) the lines bounding Le and Rp, with £' = 
{£x — l,£y — 1) intersect below the line containing ze- We then obtain E{£) C Rp n Bp f\Bi = %. In 
the remaining sub-case 0 G [vr, 57 r/ 4 ] the regions A(S'(^, u)) and A{Be) are disjoint. It follows that 
E {£) is empty. 

case 0 G (57r/4,37r/2). Point £ is now the point with the maximum y-coordinate. It then follows that 
all labels of S{£, u) lie in Bp, where £' = {£x, £y + 1). Hence, we also get E{£) c Bp. The labels of 
S{£,u) do not intersect ze, hence they are contained in LeUTe. We then have E(£) C Ber\{LeUTe) = 
Be n Le. Since 0 G (37r/4, 37r/2) it now follows that the height and width are both at most one. 
case 0 = 37r/2. All labels from S{£, u) lie in Le, all labels from Be = Re lie in Re. Hence, E{1) = 0. 
case 0 G (37r/2, 77r/4). It is again easy to show that E{£) C Re Fi Lr, where r = {£x + 1, £y + 1), and 
thus has width (at most) one. All labels from points in S{£, u) have width and height one, and are thus 
contained in Lr. Furthermore, they do not intersect Zi, from which we obtain that they are contained in 


e(0) < < 


1 //■0 = 0 
2 j/0G(O,7r/4) 

1 //■0 G [7r/4,7r/2) 

0 if 9 = 7r/2 
1 /f0 G (7r/2,37r/4) 


e(0) < < 



9 








9 = 0 
1 X 1 


OeiO^ir/^) 
1 X 2 


6 G [7r/4,7r/2) 
1 X 1 


9 = -kII 



9e (7r/2,37r/4) 
1 X 1 



A. 1 

^G (57r/4,37r/2) 9 = ^t:I2 0 G {37r/2, 77r/4) 

1x1 - 3x1 



9 G [77r/4, 27r) 
2 X 1 


Fig. 5: The region £'(£) (blue) with respect to sub-problem 5(£, n) (purple), depending on the orientation 
Q of the leaders. For each orientation we give an upper bound on the size of Ei^). 
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S{i,u) L 


(a) (b) 

Fig. 6: The region E'{t) (in pink) for the cases 9 G (57r/4, 37r/2) (a) and 9 G (37r/2, 27r) (h). In both 
cases the leader 7 p of a point p G E"{1) intersects the line segment Iz and thus subdivides E{i) 
into a left region L and a right region R. 

fi U T^. From the former we get that E{i) C L^. From the latter we get that E{£) C Tp U T^/, where 
= (ix, iy — 1). Hence, we obtain E{i) C Ri D Lr Ci Lr Ci {T£i U T^/) n B£. 

Since 9 G (3'7r/2, 77r/4) the height of E{i) is determined by i' and the intersection point p between 
Bi and L^. Trigonometry now shows that the height is at most three, 
case 9 G [Tvr/d, 27r). Similar to the previous case we get a width of at most one. The height is now 
determined by and the intersection of B£ and Ry. Since 9 G [Tvr/d, 27r) the height is at most two.n 

Corollary 7. There can be at most e{9) points in E{i) labeled internally such that their labels are 
disjoint. 

Next, we turn our attention to the points in E'{j) whose leader can intersect a label of a point in 
5(f, u). A leader of a point in B£, and thus in E'{i), can intersect a label of S{i, u) only if the labels 
intersect B£. This happens only if 0 G (57r/4,37r/2) or 0 G (37r/2, 27r). See Fig. 5. In the former case 
we thus have E' {t} = B£r\T£/ f] Li, where f ^ + 1 — tan(0 — 57r/4)), and in the latter case we 

have E^i) = n Tj n T^, where z = See Fig. 6. 

We now note that if we label the set Q C E{£) internally, then all other points in E{£) are labeled 
externally. Hence, if the leaders of the remaining points (e.g. those in E{i) \ Q) intersect with labels 
of points in S{i, u), this is already captured by the configuration involving Q. Therefore, the points in 
E{£) \ Q themselves do not define new configurations. Similarly, the points in Vx do not define any new 
configurations. 

The points in E'{£) that lie outside of E{£) can still be labeled both internally or externally. Let 
E"{£) = E'{£) \ E{£) denote the region containing these points. We now observe: 

Lemma 8. Let Q be the set of points in E(£) labeled internally, and let p G E"(£) be labeled externally. 
For all points q £ Q we have: q intersects the leader jp, or if there is a label \a, a G S{£,u), that 
intersects Xq, then it also intersects the leader jp. 

Proof. It is easy to see that any leader 7 p intersects the line segment £z, with z = {£x,£y + 1) if 
9 G (57r/4, 37r/2), and z = {£x + ^i£y) if 0 G (37r/2,27r). See Fig. 6. In both cases 7 p subdivides E{£) 
into a left region L and a right region R. At any height (y-coordinate), this left region L has width at 
most one. Hence, if there is a point q £ L labeled internally, then its label intersects 7 p. Any point q £ R 
is separated from S{£, u) by 7 p. So, if there is a label a £ S{£, u) that intersects Xq, then it also intersects 

7 p. □ 

Observation 9. Let p £ V r\ E"{£) be the point closest to the slab S{£, u), and let q be any point in 
V n E"{£). If there is a point a £ S{£, u) whose label Xa intersects 7g, then Xa also intersects 7p. 
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Observation 9 give us that there is only one relevant point in namely the point p closest to 

S{1, u). Furthermore, from Lemma 8 it follows that if p exists, then there are no relevant points in E{€) 
labeled internally. Hence, p by itself determines a configuration. We can then define fhe universe U^e of 
possible configurations of as follows: 

= {(2^) 0) I C E{i) and all labels in A(X) are pairwise disjoint} U 

W,{p})\P^E"ii)} 

Let e'{6) = 1 if there is a point in E"{1) labeled externally, and e'{9) = 0 otherwise (this includes 
the case in which E'{tj = 0). Using that the labels of points in E{tj do not contain any other points, 
together with Lemma 1 we then obtain: 

Lemma 10. The number of configurations ofBi is at most 0{\hi^E\) = 


Bounding the number of configurations of T^. Analogously to the bottom influence region E{£) in 
Bi we define a top infiuence region F{u) confaining fhe poinfs from Tu whose label can infersecf a label 
of S{i, u), and a region E'{u) confaining the points whose leader can intersect a label of the points 
in S{£, u). We observe that F{u) and F'{u) are symmetric to E(£) and E'{i) by mirroring in a line 
with slope one. We thus get similar results for F and F' as those stated in Lemmas, Corollaries, and 
Observations 6-9. So, similarly we define the universe of configurations UuF of labelings of T}. We can 
then summarize our results in the following lemma: 


Lemma 11. The number of configurations of is at most 0{\hluF\) 


f{0) < 


0 

1 

0 

1 


0 

\ 


ifO = £) 
ifOe (0,7r/4) 
if 9 G [vr/4, 37r/4] 
if9£ (37r/4,7r) 
if 9 = TT 


f{0) < 


if9 e (7r,57r/4] 
if 9 G (57r/4, 3'7r/2) 
if 9 = 37r/2 
if9e (37r/2,77r/4] 
if 9 G (77r/4, 27r) 


and 



if 9 G (0, vr/4) U (37r/4, 27r) 
otherwise. 


0(6^^) + where 



3.2 Computing an Optimal Labeling 

Let 4>(^, rt, C^e^Cuf) denote the maximum number of points in S{1, u) that can be labeled internally, 
given configurations Cie = and CuF = {^u,£u)- That is, the maximum number of points in 

S{i,u) that can be labeled internally assuming that (i) the points inli C E(£) and Tu F F(u) are labeled 
internally, and (ii) the points in Eu, and the remaining points in E {£) and F (u) are labeled externally. 
An argument similar to that of Lemma 4 then gives us the following recurrence for 4>(£, u, CiE, Cue)- 


^{£,u,CeE,CuF) = max 


^4/(5(An)), 


max 

p&S{£,u), 

CpE&Q{p/,ufiiE), 

CpF ) 


4'(iipn5(£,n)) } I 

+ ^{^,P,CiE,CpE) > ? 
+ ^{P, U, CpE, Cuf) ) I 


12 








where q{p, £, u, C^e) and ?(p, £, u, Cuf) denote the universes UpE and UpF restricted to the sets compati¬ 
ble with the labeling so far, respectively. More formally, we have 


Q{p,£,u,CeE) 


/ 

(X, X) G UpE, 


X X E{j)) n {{S{£,u) n Rp) Dll), and 

(X,X) 


'{£} if ££E"{p) 


X = < 

Si if Si cE^’ip) 



0 otherwise. 

V, y 


The function is defined analogously. Similar to the previous section we define 'I'(P) as 


^{P) 


if all labels in A{P U U U {S{£, u) \ S{£, u))) are pairwise disjoint, and their 
lx*! intersection with the leaders in r('Px U {£, w} U U U {E{£)\Xi) U {F{u)\Xu)) 
is empty, 

^ —oo otherwise. 


Computing d>(£, u, Qg, We again use dynamic programming. The size of our table is now 

0{it?\Uie\\1^uf\)- To compute the value of an entry u, Cie, Cuf), we maximize over 0{n\UpE\\i^pF\ ) 

other entries. For each such entry, we need to compute the value of 'F(X’), for some set of points 
P. We can do this in 0(n log n) time using the algorithm from Section 2. In total this yields an 
0{n‘^\UiE\\b(uF\\l^pE\\^pF\ logn) time algorithm. Next, we describe how to improve this to 0{n^ logn+ 
n^\U£E\\KuF\\l^pE\\l^pF\) time by precomputing T'. 

Fix two points £ and u, and let Cie, Cuf) = '^{Rpi^S{£, u)), given configurations C^e andC^iT’. 
We use a similar approach as in Section 2. We first compute 'h' for all points p, assuming that no other 
points above or below S{£, u) interfere with S{£, u). That is, we compute all values (0,0), (0,0)). 
This takes 0(n log n) time using the same algorithm as before. For each of the remaining pairs of 
configurations {Cie, Cuf), we hnd the “rightmost” point p in S{£, u) such that the label of p conflicts 
with CdE or CuF- It then follows that ^'{p\ Cie, Cuf) = —oo for all p' ^ LpU {p}. 

Next, we describe how we can And the “rightmost” point that conflicts with C^e in O(logn) time, 
after 0(n log n) time preprocessing. We find the “rightmost” point that conflicts with CuE analogously. 
It then follows that we can compute Cie, Cuf) for all configurations and all points p ^ S{£, u) in 
0 (nlogn + \UiE\\UuF\ + {\l^eE\ + \UuF\)logn) time in total. 

The “rightmost” point that conflicts with a configuration Cie = (X, X) conflicts with the set of 
internally labeled points X, or the set of externally labeled points in X U E{£) \ X. For both these sets we 
find the “rightmost” point conflicting with it, and return the “rightmost” point of those two. 

To And the “rightmost” point q conflicting with X, we use the same procedure as in the previous 
section. We build a range tree on the corner points of A(S(£, tt)), and use a priority range query to find 
Qr G S{£,u) whose label intersects a query label Xr,r G X. We can thus find (7 in 0(|X| logn) = O(logn) 
time. 

If the labels of S{£, u) are not contained in S{£, u) then we also need to And the “rightmost” point 
whose label intersects a leader in y = r(X U E{£) \ X). To do this we need two more data structures. 
We preprocess the edges of Z = A(S(£, u)) to allow for ray shooting queries with rays of orientation 
0. Since all edges of Z are either horizontal or vertical, and all query rays have the same orientation, 
this can be done with a shear transformation and two standard one-dimensional interval or segment trees 
T' (one for the horizontal edges of Z and one for the vertical edges of Z). We can build these trees in 
0(n log n) time [ 8 ]. This allows us to find the “rightmost” point qr whose label intersects a given leader 
7 r in 0(log n) time. To find the “rightmost” point whose label intersects a leader among all leaders in Y 
we use the following approach. We query T' with the leader 7 ^ G y closest to S{£, u), and find qr (if it 
exists). We then find the first leader 7 ^ that is hit by a horizontal rightward ray starting in r (see Fig. 7), 


13 



















Fig. 7: We can find the leaders that can intersect a label from a point in S{£, u) by a series of horizontal 
ray shooting queries. 


and recursively process 7 *. For any subsequent pair of such leaders (jr, 7s) we have that point s lies 
outside of the label Ar- Since all but one of these points lie in E(£), there are only a constant number of 
such pairs. Hence, we also need only a constant number of queries in T', each of which takes 0(log n) 
time. 

The only question remaining is how to find the next leader jg given point r. To this end we maintain 
a second data structure. We use a shear transformation such that the leaders are all vertical. We then 
build a dynamic data structure V for horizontal ray shooting queries among vertical half-lines. Such 
a data structure can be build in 0(n log n) time, and allows for 0(log re) updates and queries. See for 
example [11], although much simpler solutions are possible. We can update V for the next configuration 
in O (log re) time, since only a constant number of points change from being labeled internally to 
labeled externally and vice versa. 

Hence, after 0{n log re) time preprocessing, we can compute the “rightmost” point conflicting with 
each configuration C^e in O(logre) time. The total time required to compute 'i>'{p,C£E,CuF), for all 
points p E S{i, re), is thus 0{\l(iE\\KuF\ + (n + \Uie\ + \1^uf\) log re). 

After precomputing all values of /, the dynamic programming takes 0{n^\L{iE\\£^uF\\i^pE\\£^pF\) 
time. We thus obtain a total running time of 0(re^logre + n^\U£E\\^uF\\i^pE\\£^pF\))- Using that 
\U£e\ and \UpE\ are both at most 0(re®'i^i + (Lemma 10), and that \Uuf\ and \UpF\ are at most 
0{n^ (^) -)- j/t®)) (Lemma 11) we can rewrite this to 0(re^(logre + i{n, 5, 6))), where r(re, <5, 6) is a term 
that models how much the subproblems can influence each other. We have 

i{n,5,e)= +n^e'(e)+f'(0)sm +„e'(0)+2/'(^^)^e(e) 

J^jf'(e)+ne) je(0)+/(0) _^j^2e'(e) ^ 2 /( 0 ) +^ 2 / 70 ) ^2e(0) 

je(e)+2/(0)_^^/'(0) ^2e(0)+/(0)_^ ^2e(0)-r2/(0) _ 

For 5 = 0(n), this gives us a worst case running time varying between 0(n^ log re) and 0(re^®). We 
conclude: 

Proposition 12. Given a set V of n points and an angle 9, we can compute a labeling of V that 
maximizes the number of internal labels in 0(re^(logre + i{n, 5,6))) time and 0{inf i{n, 6, 9)) space, 
where 6 = min{re, l/dmin}/or the minimum distance dmin in V, and i{n,8,9) models how much 
subproblems can influence each other. 

3.3 An Improved Bound on the Number of Configurations 

The analysis above, together with the fact that e{9) and f{9) are both at most three, gives us a worst case 
running time of 0{n^^). We now study the situation a bit more carefully, and show that the number of 


14 






Fig. 8: S{£, u) is incompatible with any point a £ Cl Lr HTr (the orange region). 

interesting configurations is much smaller. More specifically, that we can replace e{6) and Sf{0) by a 
quantities e*{9) and f*{9) that are both at most one. This significantly improves the running time of our 
algorithm. 

We start by observing that for some points q in E{£), the subproblem u) does not have a labeling 
compatible with q, irrespective of whether q is labeled internally or externally. Let Q{tj denote the set 
of such points. It follows that we can restrict ourselves to labelings, and thus configurations, that do 
not contain points from Q{1). Let = {{X,E) \ {X,E) £ IAie A X C !"(£)} denote this subset of 
configurations, where y{€) = E{i) \ Q{£). 

Lemma 13. For every configuration (X, X) £ ^ most two. 

Proof. By Corollary 7 there are at most e{9) points in E{i) that can be labeled internally simultaneously. 
Hence |X| < e{£). For 9 £ [0,37r/2] U [77r/8,27r] we have e{9) < 2, and thus the lemma follows 
immediately. For 9 £ (Ttt/S, 27r) we have e{9) < 3. We prove this remaining case by contradiction. 

Assume that 9 £ (77r/8, 27r), and that X = {a, b, c}. Since the region E{i) has size 1 x 3, it follows 
that one of these points, say point a, lies in the triangular region B^n LrC] Tr, with r = [ix + + 1). 

See Fig. 8 . Let p be a point in S{1, u) whose label intersects Aa- Using that 9 £ (Fvr/S, 27r) it follows 
that p is to the bottom-left of a. Since the labels are to the top-right of a point, we then obtain that a £ Xp. 
Therefore, the leader 7 ^ also intersects Xp. So, both the label and the leader of a interfere with S{£, u), 
and thus a ^ X{£) and thus also not in X C X(^). Contradiction. □ 

Lemma 14. Let p be a point in S{£, u), and let (X, E) £ IL^e- label ofp intersects at most one label 
Xq of a point q £ I. 

Proof. When |X| < 1, the lemma is trivially true. By Lemma 13, we otherwise have |X| < 2 and 
9 £ (0,7r/4) U (37r/2,27r). Let X = {a, 6 }. For the case 9 £ (0,7r/4), assume w.l.o.g. that b is the 
leftmost point. Since a and b are both in X, their labels are disjoint, and thus we have hx + '£ < ax (See 
Fig. 9(a)). Furthermore, we have Py > £y > by. Since Xp intersects Xj,, butp 0 Xi, this means px < bx. 
Since Xp has width one, it thus cannot intersect Aa- 

For the case 9 £ (37r/2, 27r) we use a similar argument. Assume w.l.o.g. that a is the topmost point, 
and thus Oy > by + 1. See Fig. 9(b). Since p £ S{£, u), Xp intersects Aa, and a ^ XpWe have that p is to 
the top-left of a. And thus, py > Oy > by + 1. The label of b has height at most one, and thus cannot 
intersect Xp. □ 
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(a) 


(b) 




Fig. 9: Point p £ S{£, u) can intersect only one label of a point in E{1). 


Lemma 15. Let C = ({a, 6}, £) G where a is to the right ofb ifO G (0,7r/4), and to the top ofb if 
9 G (37r/4, 2ti). Point b uniquely determines a. 

Proof. We prove this by contradiction. Assume that there is another configuration C = ({o', b},£') G 
such that a' £ Y(£) is to the right of b in case 9 £ (0, 7 r/ 4 ), or above b in case 9 £ ( 37 r/ 4 , 27r). 

The points a and b can influence the labeling of S{£, u). Hence, their labels intersect with labels of 
points in S{£, u). By Lemma 14 Xa and Af, cannot both intersect the same label of a point in S{£, u). So, 
let p be the point whose label intersects Aq, and let q be the point whose label intersects Xb- See Fig. 10. 

We start with the case 9 £ (0, 7 r/ 4 ). Point a is to the right of b, and Aa and A;, are disjoint. Hence, 
(lx > bx + 1. Since Ap intersects Xa, and p 0 Xb, it follows that p lies above Xp. Hence py > by + 1. 
Again using that Xp intersects Xa, it then also follows that Oy > by. Using the same argument we get 
a'y > by. Now in the configuration C, point o' is labeled externally. However, this means its leader 
intersects A;,. Hence, C 0 Contradiction. 

For the case 9 £ ( 37 r/ 2 , 27r) we have that 6 ^ + 1 < Oy. As in the proof of Lemma 13 we have that 
Oy < £y + 1, and hence by < £y. Using that Xq does not intersect 7 ^, we have qy > £y. Since q ^ Xbit 
follows that qx < bx. In configuration C, Xa' and Xb are also pairwise disjoint. It then follows that o' is 
to the top-right of b. In the configuration C = {{a,b},£), point o' is labeled externally. However, this 
means its leader intersects Xb. Hence, C 0 Contradiction. □ 

Corollary 16. The number of configurations in is at most where e*{9) = 

min{l, e(0)}. 

We can again use a symmetric argument for the number of configurations in Tu. Figure 11 gives a 
graphical summary of these results. We now redefine r(n, <5, 9) as 

tin, 6,9)= +^2e'(e)+/'(e)^r W +„e'(e)+2/'W^e*(e) 

j^ne'(e)+r(e) ^e*(e)-r/*(0) ^2/*(0) ^2e*(e) 

+„e'(0) §e*ie)+2r{e)_^^fi9) ^2e*(0)+/*(0)_^ ^2e*(0)-r2/*(0)_ 

and obtain the following result, which is at most 0{'nJ) for 6 = 0{n): 
























Fig. 10 : The points a, a', and b in E{1) and the label(s) in S{1, u) they intersect. (Fig (a) is not on scale.) 



37r/2 

Fig. 11: A depiction of the upper bounds on e*,e', s*F, and /' as a function of 9. 


Theorem 17. Given a set V of n points and an angle 6, we can compute a labeling ofV that maximizes 
the number of internal labels in 0{n^{logn + i{n,5,6))) time and 0{'n?i{n, 6,9)) space, where 
6 = min{n, l/dminl/or the minimum distance dmin in E, and L{n, 6, 9) models how much subproblems 
can influence each other. More formally, 

i{n,5,9) = +n^e'{e)+ne)^r(e) ^^e'{e)+ 2 ne)^e^ie) 

(0) (e)+2/* (0) ^26* (0)+/* (e) _p j2e* (0)+2/* (0) ^ 

where e*{9), f*{9), e' {9), and f'{9) are all at most one. 

4 Extensions 

So far, we have considered a stylized version of the question we set out to solve. In this section we discuss 
how our solution may be adapted and extended, depending on the exact requirements of the application. 

4.1 Optimizing the direction 

Rather than fixing the direction for the leaders in advance, we may be willing to let the algorithm specify 
the optimal orientation that maximizes the number of points that can be labeled internally. Or, perhaps 
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(b) 6 — 7r/3 


(C) 6» = 77r/4 


(a) e = 0 


(d) e = ll7r/6 


Fig. 12: Examples for routing the outer leaders and placing the external labels for different slopes. 


we wish to compute a chart that plots the maximum number of internally labeled points as a function of 
the leader orientation 9, leaving the final decision to the judgement of the user. 

In both scenarios, we need to efficiently iterate over all possible orientations. We adapt our method 
straightforwardly. Let Q be the set of all 4n corner points of all potential labels. For every pair p,q ^ Q 
consider the slope 9p^q of the line through p and q. All values 9p^q partition all possible angles into 0{'n?‘) 
intervals. For all values 9 in the same interval J, any leader 7 ^ intersects the same set of potential labels, 
so the optimal set of internal labels is constant throughout J. We compute it separately for each interval. 

By applying Theorem 17, we achieve a total of 0(n^ • n^(logn + i{n,5,9))) = 0(n^(logn + 
i{n, 6■,9))) time to compute the optimal labelings for all orientations, or to optimize the orientation by 
performing a simple linear scan. 

4.2 Routing the outer leaders 

Once the core combinatorial problem of deciding which points have to be labeled internally is solved, 
it remains to route the outer leaders and place the external labels. Since our goal in this paper is to 
maximize the number of internally labeled points, we are only interested in finding a valid labeling, in 
which neither labels nor leaders intersect each other. Let us assume that 9 G [0, 7 r/ 2 ] U [ 37 r/ 2 , 27r], i.e., all 
external labels are oriented to the left. The case of labels oriented to the right is symmetric. We consider 
the leaders in counterclockwise order around the boundary of M. and place them one by one starting with 
the topmost leader. The first label is placed with its lower right corner anchored at the endpoint of its 
inner leader. For all subsequent labels we test if the label anchored at the endpoint of the inner leader 
intersects the previously placed label. If there is no intersection, we use that label position. Otherwise, 
we draw an outer leader extending horizontally to the left starting from the endpoint of the inner leader 
until the label can be placed without overlap. Obviously this algorithm takes only linear time. Figure 12 
shows the resulting labelings for four different slopes. We note that depending on the slope 9 of the inner 
leaders other methods for routing the outer leaders might yield more pleasing external labelings. This is, 
however, beyond the scope of this paper. 

4.3 Non-square labels 

Square labels are not very realistic in most map-labeling applications. Their use is justified by the 
observation that if all labels are homothetic rectangles, we can scale the plane in one dimension to obtain 
square labels without otherwise changing the problem. Nonetheless, reality is not quite that simple, for 
two reasons: firstly, the scaling does alter inter-point distances, so if we wish to parametrize our solution 
by dmin we need to take this into account. Secondly, in real-world applications, labels may arguably have 
the same height, but not usually the same width. 

If all labels are homothetic rectangles with a height of 1 and a width of w, scaling the plane by a 
factor 1/win the horizontal direction potentially decreases the closest interpoint distance by the same 
factor. Now, the number of points in a unit-area region that do not contain each other’s potential labels 
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is bounded by dw, immediately yielding a result of 0(n^(log n + i(n, 6w, 6))) using exactly the same 
approach. 

When all labels have equal heights but may have arbitrary widths, we conjecture that a variation of 
our approach will still work, but a careful analysis of the intricacies involved is required. If the labels 
may also have arhitrary heights the problem is open. It is unclear if there is a polynomial time solution in 
this case. 

4.4 Obstacles 

In this paper, we have considered only abstract point sets to be labeled, using leaders that are allowed to go 
anywhere, as long as they do not intersect any internal labels. While this is justified in some applications 
(e.g., in anatomical drawings, it is common practice to ignore the drawing when placing the leaders, as 
they are very thin and do not occlude any part of the drawing), in others this may be undesirable (in 
certain map styles, leaders may be confused for region boundaries or linear features). As a solution, we 
may identify a set of polygonal obstacles in the map, that cannot be intersected by leaders or internal 
labels. 

In this setting, obviously not every input has a valid labeling: a point that lies inside an obstacle can 
never be labeled, or obstacles may surround points or force points into impossible configurations in more 
complex ways. Nonetheless, we can test whether an input has a valid labeling and if so, compute the 
labeling that maximizes the number of internal labels in polynomial time with our approach. 

The main idea is to preprocess the input points in a similar way as in the beginning of Section 3. 
Whenever a point has a potential leader that intersects an obstacle, it must be labeled internally; similarly, 
whenever a point has a potential internal label that intersects an obstacle, it must be labeled externally. If 
we include such “forced” leaders or labels into our set of obstacles and apply this approach recursively, 
we will either find a contradiction or be left with a set of points whose potential leaders and potential 
internal labels do not intersect any obstacle, and we can apply our existing algorithm on this point set. 

The same approach may be used for point sets that are not in general position: if we disallow leaders 
that pass through other points, they are forced to be labeled internally. Note that this again may result in 
situations where no valid labeling exists. 
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